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The interaction of spatially extended heavy hadrons is investigated in the framework 
of lattice QCD with dynamical quarks. In addition to the baryon-baryon potential results 
for the baryon-antibaryon and for the meson-meson system are presented. It is shown 
that the expected dipole forces have a very short range and that sea quarks play a minor 
important role. 
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QCD as a quantum field theory has due to the uncertainty principle a creation pro- 
bability of virtual gluons and quarks. As a consequence the nucleon-nucleon forces are 
mediated by gluon exchange between the constituent quarks for short distances whereas 
for longer distances the production of quark-antiquark pairs is the dominating mechanism 
[1]. The continuous advance of the techniques of lattice QCD makes the study of systems 
consisting of a few valence quarks in the presence of sea quarks feasible [2]. This paper 
investigates the influence of both gluons and dynamical quarks on static hadron-hadron 
systems. The sea quarks constitute the meson exchange as expected by the Yukawa 
theory. 

In our approximation of QCD the valence quarks of a hadron are restricted to ful- 
fill the static Dirac equation. The gluons are treated as SU(3)-Maxwell fields U^^ and 
the sea quarks as Dirac fields t/j^^t/j^ in Kogut-Susskind discretization. The free energy 
F(ri, . . . ,rjv) of a system of N quarks and antiquarks in a gluonic and fermionic field 
is defined as a thermodynamical expectation value and expressed by the Feynman path- 
integral for the product of N Polyakov loops L[fi) [3] 



S{U^t/j^t/j) = S'G(t^) + 5'i?([/, ?/',?/') is the total action of the system with temperature r. For 
the gluonic action Sg we use Wilson's plaquette action with the inverse gluon coupling 
/3 = 6/(/^ and for the fermionic action Sf we employ the Kogut-Susskind formulation 
with nj fiavors. The time evolution of the valence quarks is taken as the static quark 
propagator, the so-called Polyakov loop L[f) = 1/3 Tr Hil:'! ^^^4=4(^5 ^O- construct a 
quark wave function ^(r) by distributing the static quark charges in Gaussians 



over a sphere with radius R on the lattice. Since the results turn out to depend only 
weakly on the width a we set (2) to a uniform distribution. For the description of 



exp ( - 7f;F{ri, ■ ■ • ,nv)) 



fV[U]V[tlj,ij]e-^\-U'^'^^ 
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baryons we take the product of three quark wave functions and for mesons correlated 
wave functions of quark and antiquark. Because the static quarks carry no spin and flavor 
all baryons and mesons are degenerate. The total wave function of the hadron-hadron 
system is constructed to be a product of Gauss functions with two centers separated by 
some distance d. 

To compute the path integral (f) we perform fOOOO Monte Carlo iterations for QCD 
without dynamical fermions and f 000 iterations in the presence of fermions. We take a 
space-time lattice of size 8 X 8 X f6 X 4 with periodic boundary conditions. In the pure 
gluonic case we choose /3 = 5.6 which corresponds via the renormalization group equation 
to a lattice constant of a 0.25 fm and to a spatial extension in z direction of about 
2 fm. For full QCD at /3 = 5.2 (/3 = 4.9) we set the number of flavors nj = 3 (nj = 4) and 
the dynamical quark mass m = 0.1. This choice leads to comparable lattice constants all 
lying in the confinement regime. 

We calculate the total energy of the hadron-hadron system according to (1) and sub- 
tract twice the energy of a single hadron. In Fig. la we show the relative potential energy 
of the two-baryon system with baryons of radius R = 1 for varying distances d between 
the two clusters. We find an attractive potential with an effective range of about 2R. The 
reason for the attractive forces between the two baryons is the Coulomb plus linear type 
of potential between quarks which tries to attract the two three-quark clusters. When 
the two baryons become separated their colors are saturated and they evolve to isolated 
color singlets which do not interact. At least on the lattice no long range forces can be 
resolved. Taking dynamical fermions into account we obtain no drastic effect. This is also 
observed in mass spectrum calculations where sea quarks give rise to only a 10 per cent 
effect. In Fig. lb we display the two-baryon system with baryon radius R = 2. We find 
attractive interactions only in the overlap region and practically no fiavor dependence. 

We now turn to the baryon-antibaryon system in Fig. Ic and again see an attractive 
interaction in the overlap regime. The difference between the baryon-baryon and baryon- 
antibaryon potentials might be interpreted qualitatively by the fact that the baryon- 
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antibaryon system can build three mesons which is energetically favored. Switching on 
dynamical fermions we observe practically no effect. Increasing the radius of the hadron 
sphere to i? = 2 we recognize in Fig. Id more or less the same behavior as for R = 1. 
The static meson-meson system is depicted in Figs. I e-f for radii R = 1 and i? = 2, 
respectively, and behaves similar to the previous systems. 

To conclude, this study deals with a truncated Dirac equation restricting the valence 
quarks to static color charges. The construction of the quark orbitals leads to baryons 
of reasonable size but takes no kinetic energy into account. It turns out that the gluon 
exchange dominates and sea quark effects are small. The aim would be to evaluate the 
S matrix for scattering of two nucleons consisting of three quarks each. This requires 
the calculation of a four-point hadron- Green-function being equivalent to a correlation 
function of six quark propagators in the QCD path-integral. There has been some pro- 
gress in computing meson scattering-amplitudes from the finite volume dependence of the 
bound-state energies [4]. A method to extract effective potentials from the two-hadron 
spectrum was proposed recently [5]. 
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Fig. 1: Baryon-baryon potentials (a,b) as a function of center of mass distance d for pure 
gluon exchange and with sea quarks for baryons with radii R = 1 and i? = 2 in lattice 
units. The interaction is attractive and takes place mainly in the overlap region. The 
baryon-antibaryon potentials (c,d) turn out to be deeper and the meson-meson potentials 
(e,f) are shallower. Dynamical quarks have no pronounced effect. Error bars are in the 
size of the symbols. 
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